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Abstract. We study transport through a double quantum dot with interdot hopping t, 
intradot repulsion U and interdot repulsion 17', using the numerical renormalization group 
(NRG) method. At half-filling, the conductances in two-terminal series and four-terminal 
parallel configuration are calculated via two phase shifts for quasi-particles of double quantum 
dot connected to two noninteracting leads with hybridization strength F. For small values of 
t/T and U' /U, conductance in the two-terminal series configuration is suppressed to almost 
zero. In this region, plateau of conductance in the four-terminal parallel configuration appears 
and almost reaches a unitary limit value Ae^/h of two conducting modes. For large values of 
t/r or U' /U, both conductances are suppressed to almost zero. The conductance in the two- 
terminal series configuration almost reaches a unitary limit value 2e^//i only around cross-over 
regions of electron-configuration in double quantum dot. Through the behavior of the local 
charge and some thermodynamic quantities, we discuss the relation between transport and 
electron-configuration. 



A half- filled double quantum dot with interdot repulsion has been theoretically studied [HEIE] 
and some interesting phenomena, for example, quantum phase transition at a critical interdot 
repulsion [1], have been predicted. In this paper, focusing on two kind of conductance through a 
half-filled double quantum dot with interdot repulsion, we discuss the relation between transport 
and electron-configuration of double quantum dot through the behavior of the local charge 
and some thermodynamic quantities, on the basis of numerical renormalization group (NRG) 
calculation. 

We consider two-site Hubbard model with interdot repulsion, which is connected to two non- 
interacting leads at the left(L) and right(R) by the symmetrical tunneling matrix elements v, as 
illustrated in Fig. [T] (a). 

The Hamiltonian is given hy H = + Hj^m + -ffiead with 
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Figure 1. Schematic picture of (a) series configuration with two terminals and (b) parallel 
configuration with four terminals 



where dif, annihilates an electron with spin a at site i in the double quantum dot, which 
is characterized by the interdot hopping matrix element t, onsite energy e^, intradot repulsion 
\J and interdot repulsion U' . In the lead at v (= L, i?), the operator c^^^^ creates an electron 
with energy corresponding to an one-particle state (pkuif)- The linear combinations of the 
conduction electrons ipj^^ and V/jo- mixed with the electrons in dots labeled by z = 1 and 2, 
respectively, where tp^^^ = ^k'^kua^ku{'''y) and r,j is the position at the interface in the lead 
V. We assume that the hybridization strength F = 'KV^Ylfk \(t>ku{'rv)\^ ^{^^ — ^^u) ^ constant 
independent of the frequency uj and i^, and take the Fermi energy ^ to be /u = 0. 

Our system has an inversion symmetry, so we can introduce the even and odd-parity orbitals 
as follows; 

Two phase shifts Se and 5o of the quasi-particles with even and odd parities characterize a local 
Fermi-liquid behavior of the whole system described by H at low temperature. We can obtain 
5e and 6o from the fixed-point eigen values of NRG. From two phase shifts defined with respect 
to the system described by H, we can deduce not only the conductance gg in the two-terminal 
series configuration illustrated in Fig. [T] (a) but also the conductance gp in four-terminal parallel 
configuration illustrated in Fig. [1] (b) at T = as follows; 
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We calculate series conductance gs, parallel conductance gp, electron number rig = Yjai^t'^o-) 
in even-parity orbital defined by Eq.Q and some thermodynamic quantities at half-filling 
(e^ = -U/2 - U') for 10"^ < t/T < 10 and < U'/U < 1.5. In this paper, we fix U/T to 
he U/T = 15. In Fig. [2l (a) gs, (b) gp and (c) rig are plotted as functions of t/T and U'/U. First 
of all, to make discussion clear, we consider three limit case; case 1: t/T and U'/U 
, case 2: t/T oo and U'/U , and case 3: t/T and U'/U oo. In Fig. [2]^a), we 
illustrate schematic picture of electron-configuration in double quantum dot for these three limit 
case. In the case 1, because of no interaction between two dots, our system is decoupled to two 
same systems composed of a single dot connected to a single lead. At half-filling, each single 
dot is occupied by one electron to avoid intradot repulsion U, as shown in the picture [case 1] in 
Fig. [2l||a). Schematic picture [case 2] in Fig. [2l||a) represents that two electrons occupy the even- 
parity orbital stabilized by large interdot hopping t in the case 2. So, we can expect that the 
electronic states of the ne — 2 region in Fig. E^c) are similar to the state of the case 2. The case 
3 have been studied in detail by Galpin et al. using NRG|1, Jj. Either of two dot is occupied by 



two electrons to avoid strong interdot repulsion U' which is much larger than intradot repulsion 
[/, as shown in the picture [case 3] in Fig. [2ja). Therefore, a degenerate configuration of two 
electrons results in non- Fermi liquid ground state with ln2 residual entropy. From Fig. [2]^a) and 
above discussion, we find that the series conductance almost reaches a unitary limit value 
le^ jh only around cross-over regions of electron-configuration in double quantum dot. Next, 
we discuss the parallel conductance g-^ in three limit cases. From the similar discussion for 
two-terminal series configuration, we can easily understand the schematic pictures of electron- 
configuration in four-terminal parallel configuration for three limit case, shown in Fig. [2][b). For 
the case 1, our system is decoupled to two same systems composed of single-dot connected to 
two leads. In each system, one electron occupies single dot as shown in the picture [case 1] in 
Fig. mb), and causes Kondo effect. Therefore, the parallel conductance reaches a unitary 
limit value equal to ^^e^ jh of two conducting modes. From Fig. E^b), we indeed find that there 
is plateau of the parallel conductance (7p, which almost reaches a unitary limit value 4e^//i, for 
small values of t/F and U' /U . The even-parity bonding orbital of double-dot is fully occupied 
by the two electrons in the case 2. In the case 3, two electrons fully occupy either of two dots 
and no electron occupies the other, as shown in the picture [case 3] in Fig. [2)||b). Therefore, the 
parallel conductance in both case 2 and 3 is 0. This corresponds to the fact that the parallel 
conductance for large values of t/T or U' /U is suppressed to almost zero, as shown in Fig. [2]J^b). 
Fig.[2^a) is consistent with the phase diagram obtained by Mravlje et al. using Gunnarsson and 
Schonhammer projection-operator method [3j, apart from details. 

We consider also impurity entropy S and impurity spin susceptibility x calculated using NRG 
eigen-energy. In Fig. [3l we show the temperature dependence of (a) impurity entropy S and 
(b) xT for three points in t/T - U' /U plane shown in Fig. El^c); {t/T,U'/U) = (0.0001,0.2), 
{t/T,U'/U) = (10,0.2) and {t/T,U'/U) = (0.0001,1.5). For {t/T,U'/U) = (0.0001,0.2), 
impurity entropy shows steps of ln4 at T/D ~ 10~^ and vanishes as temperature decreases. 
There is a corresponding shoulder in impurity susceptibility. It implies that the ln4 step 
originates from the degree of spin freedom. We can understand these thermodynamic behavior 
from consideration of the case 1. Each 1/2 - spin in single dot of two same systems for the 
case 1 is screened by Kondo effect in the same manner. Therefore total entropy 2xln2=ln4 
decays to all at once. The impurity entropy and susceptibility for {t/T,U' /U) = (10,0.2) 
decrease and vanish at high-temperature T/D ~ 10^^. This is because, in the case 2 assigned 
for (t/F, [/'/[/) = (10,0.2), the even-parity bonding orbital of double-dot stabilized by large 
hopping t is fully occupied by the two electrons and the freedom of spin and orbital are quenched 
at T/D ~ 10"^ For {t/T,U'/U) = (0.0001,1.5), a plateau of ln2 is seen in the impurity 
entropy as temperature decreases. The impurity spin susceptibility dose not change around 
temperature region 10~^^ <T/D< 10~^^ where the impurity entropy begins to drop from ln2 
to 0. Therefore, we find that the ln2 entropy is made up of the degree of orbital freedom. The 
ln2 plateau corresponds to non- Fermi liquid ground state with residual entropy in the case 3. 
But for {t/T,U' /U) = (0.0001, 1.5), a finite t results in Fermi liquid ground state and thereby 
the impurity entropy vanishes as temperature is lowered. 

In summary, we have studied transport through a half-filled double quantum-dot with interdot 
repulsion in two-terminal series and four-terminal parallel configuration, the electron number 
in even-parity orbital, impurity entropy and spin susceptibility, using NRG method. Through 
the behavior of these quantities, we have presented the relation between transport and electron- 
configuration of double quantum dot in wide parameter region. 
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Figure 2. (a) series conductance Qs/ {2e^ /h)^ (b) 
parallel conductance gp/{2e'^/h) and (c) electron 
number in even-parity orbital plotted as functions 
t/T and U' /U. In (a) and (b), schematic picture of 
electron-configuration for case 1, case 2 and case 3 
are presented. 
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Figure 3. The temperature dependence 
of (a) impurity entropy S and (b) 
xT for three points in t/T — U'/U 
plane shown in Fig. ^c);{t/T ,U' /U) = 
(0.0001,0.2), (t/r, [/'/[/) = (10,0.2) and 
[t/T, U'/U) = (0.0001,1.5). 
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